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ABSTRACT 

Modified  grey-scale  image  objects  classification  on  the  basis  of  variance  analysis  of  their  vectorial  models 
is  suggested.  After  calculating  the  generalized ,  the  intergroup  and  the  average  of  particular  dispersions 
the  variance  ratio  of  vectorial  models  of  analyzing  objects  is  computed  and  classification  on  the  basis  of 
threshold  values  of  their  similarity  is  made.  Experimental  classification  data  are  cited. 


INTRODUCTION 


Similarity  measure  is  one  of  the  important  instruments  of  data  analysis  in  matching,  classification, 
similarity  evaluation  and  correlation  of  analyzing  objects  [1].  Wide  application  of  the  correlation 
coefficient  as  a  measure  of  similarity  of  analyzing  objects  becomes  above  all  clear  due  to  its  stability  as 
an  image  recognition  characteristic  both  under  gross  and  small  volume  of  sampling  [2].  It  should  be 
stated  that  the  correlation  coefficient  is  near  true  as  few  as  there  is  a  gross  number  of  independent 
observations  subject  to  the  same  near-normal  distribution.  In  cases  with  limited  number  of  experimental 
data  and  their  non-linearity  it  is  recommended  [2]  to  use  correlation  ratio  (CR)  -  tj. 


?t  = 


(i) 


—  2  2 

where  8 1  -  intergroup  (intragroup)  dispersion;  G  -  generalized  dispersion. 


1.  CORRELATION  RATIO  MODIFICATION 

Weights  of  vectors  of  the  tracing  contour  described  similarly  the  chain  code  (table  1)  would  be  used  as 
invariant  information  features. 

Table  1 


Normalized  vectorial  models  of  analyzing  objects 


Object  number 

Vector  weights 

Rotation  angle,  a 

1 

2, 2, 1,2, 4, 5, 6, 5, 5, 5, 7, 8, 8  (13) 

0° 

2 

2, 1,2, 2, 4, 5, 6, 5, 5, 5, 7, 8, 8  (13) 

0° 

3 

2, 2, 1,2, 4, 5, 5, 6, 5, 5, 7, 8, 8  (13) 

0° 

Weights  of  vectors  of  the  object  #1  used  as  etalon  features  are  labeled  as  A1 .  Features  of  analyzing 
objects  #2  and  #3  (table  1)  are  correspondingly  labeled  as  B2[i ]  and  B3[i].  Visualization  of  3  object 
contours  and  superimposed  image  of  these  synthesized  objects  are  shown  on  fig  3. 
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Fig.l.  Sub-superimposed  (a)  and  superimposed  (b)  image  of  3  synthesized  objects 
The  average  of  particular  dispersions  [3]  for  objects  #2  and  #3  is  correspondingly 


a2A(B2)  =  =2l  =  0,15384615; 

n  13 

.  Y(4']-£3m)2  2 

a2A(B3)  =  ^  LJ - =  —  =  0,15384615. 

n  13 


Generalized  dispersion 


2  12  342  60  *-r>n<ai'7i£ 

<7 ,  =  — - A,  = - =  5,00591716. 

n  13  13 


Intergroup  dispersion  for  objects  #2  and  #3 

72  -2  =2  =4,85207101;  5 2  = 


2  —2 
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=  4,85207101. 


Correlation  ratio 


1 4,85207101 
v  5,00591716 

1 4,85207101 
v  5,00591716 


0,98451366; 

0,98451366. 


Calculated  CR  of  analyzing  objects  B2  and  53  relative  to  CR  of  the  etalon  are  equal  although 

descriptions  of  their  contour  representation  are  different  (table  1).  Analysis  of  vector  weights  and 
experiments  to  qualify  correlation  ratio  have  shown  [5,  6]  that  it  is  necessary  to  inset  ranking  of 
differential  values  (A[i]  -  Bl  \f[)k  .  Insertion  of  k-coefficient  for  differential  values  ranking  of  vectors  is 

the  first  modification  of  the  utilized  correlation  ratio  [6,  9].  Differential  values  ranking  could  be  made 
using  the  following  mathematical  equations: 


KAlQ  +  BJi]) 

V  n(s>  +  si ) 


(2);  k2  = 


4] !  bM 
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(3); 
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where  i  -  running  number  of  the  analyzing  object  vector  whose  differential  weight  is  other  than 
zero; 

n  -  number  of  vectors  of  etalon  or  the  analyzing  object  (the  most  of  them  is  chosen); 

5j  -  sum  of  weights  of  vectors  of  the  etalon; 

S[  -  sum  of  weights  of  vectors  of  analyzing  object  (7). 

Second  modification  of  CR  is  bound  up  with  the  requirement  fulfillment  about  the  similarity  measure 
symmetry  [2]  7(AB)=7(ba )•  Modified  CR  (MCR)  meeting  mathematical  requirements  of  similarity 
measure  [2]  is  given  by  the  following  [5]: 


*1m  = 


(5) 


where  cr^  =  o\  +  <j2b  -  sum  of  generalized  dispersion  of  the  etalon  (A)  and  the  analyzing  object  (B/); 

=  Sj  +  SB  -  sum  of  intergroup  dispersion  of  the  etalon  (A)  and  the  analyzing  object  (B /). 

Experiments  gave  preferences  to  the  ranking  coefficient  selection  and  the  add-on  of  lacking  direction 
vectors  up  to  the  equinumerous  set  of  features  by  mathematical  expectation  meanings  of  their  vectorial 
models: 


rl(AiB2)  ~rl{ b2,a)  —0,99437419;  ~71  (b3,a)  0,99209212. 


2.  VARIANCE  RATIO  FORMATION  OF  VECTORIAL  MODELS  OF 
ANALYZING  OBJECTS. 

On  a  basis  of  MCR  utilized  as  a  similarity  measure  of  analyzing  objects  boundary  representation  on  grey¬ 
scale  images,  we’ll  consider  the  following  variance  ratios: 

•  intergroup  and  generalized  dispersions  (5)  -  correlation  ratio; 

•  generalized  dispersion  to  the  sum  of  the  generalized  and  the  average  of  particular 
dispersions 


'll  = 


2  — 2 

°R,  +  a  I 


(6) 


•  intergroup  dispersion  to  the  sum  of  the  generalized  and  the  average  of  particular 
dispersions 


73  = 


2  — 2 

°Rt  +  °l 


(7) 


aimed  for  evaluation  of  a  similarity  measure  and  subsequent  objects  classification. 


3.  EXPERIMENTAL  VERIFICATION  OF  AN  EQUINUMEROUS  SET  OF 
FEATURES  FOR  OBJECTS  CLASSIFICATION. 

To  process  a  grey-scale  remote  sensing  space  image  on  a  computer  it  is  necessary  to  digitize  it.  A  digital 
image  of  certain  locality  is  depicted  on  the  figure  2.  Optical  density  (image  density)  value  and  its  outside 
configuration  (geometry)  characterize  a  grey-scale  image  of  the  locality.  A  contour  carries  the  maximal 


RTO-MP-SET-094 


10-3 


UNCLASSIFIED/UNLIMITED 


UNCLASSIFIED/UNLIMITED 


Possible  Procedure  Modification  and  a  System 
Likeness  for  Object  Identifying  on  Remote  Sensing  Images 


ORGANIZATION 


information  about  an  object  image  and  significantly  enables  to  decrease  the  computer  input  of  information 
content  during  selective  inputting  but  particularly  affected  to  brightness  and  geometric  distortions.  In  this 
context  utilization  of  optimum  intervals  of  brightness  threshold  values  and  local  anisotropic  filtration  of 
search  zones  of  the  nearest  point  enables  to  get  more  higher  boundary  representation  of  grey-scale  image 
objects  and  forms  its  vectorial  model  in  a  boundary  tracking  process. 

The  digital  image  after  its  filtration  and  boundary  extraction  utilizing  the  concurrent  activities  in 
displacement  and  summation  during  the  image  convolution  with  mask  transform  [5]  is  depicted  on  the 
Figure  3.  In  process  of  boundary  extraction  the  vectorial  models  of  boundary  representation  of  analyzing 
object  are  formed. 


Fig.  2  Digital  image  of  the  locality  Fig.  3  Digital  image  processing  effect 

To  increase  accuracy  and  trustworthiness  of  contour’s  readout  on  grey-scale  object  images  in  the  United 
Institute  of  Informatics  Problems,  Belarus  it  was  developed  a  system  providing  automatic  selective 
readout  and  boundary  representation  [8].  The  subsystem  for  boundary  representation  (Fig.  4)  contains 
Scanner  (S),  Photoelectric  converter  (PhC),  Analog-digital  conversion  unit  (ADCU),  Memory  unit  (MU), 
Boundary  points  selection  unit  (BPSU),  Control  unit  (CU)  Searching  scan  unit  (SCU),  Graphics  tracking 
unit  (GTU),  Threshold  value  formation  unit  (TVFU)  and  Commutator  (C) 


Fig. 4  Block  diagram  of  object  edge  recognition  subsystem  on  remote  sensing  images 
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Normalization  of  analyzing  objects  vectorial  models  by  scale  and  orientation  [6]  enables  to  form  invariant 
information  features.  Normalized  vectorial  models  of  7  analyzing  objects,  first  of  which  was  set  for  the 
etalon  are  located  in  table  2. 

Adding  the  lacking  vectors  weights  up  to  the  equinumerous  set  of  features  by  mathematical  expectation 
meanings  and  using  the  ranking  coefficient  we  compute  variance  ratios  (5)  and  (6)  for  the  locality 
situated  analyzing  objects  by  the  worked  out  software  module  TK07,  and  classify  them  by  threshold 
values  of  similarity  measure  [7]. 


Table  2 

Normalized  vectorial  model  of  initial  image  objects 


No 

Vectors  weights  {number  of  vectors} 

Rotation  angle, 
a 

1 

1,1, 1,1, 3, 3, 3, 3, 3, 5, 3, 3, 3, 5, 5, 5, 5, 7, 5, 7, 7, 7, 1,7, 7, 7, 1,7  {28} 

0° 

2 

1,1, 1,3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 5, 5, 5, 5, 5, 7, 7, 7, 7, 7, 7, 7, 7, 1,1, 7, 7  {30} 

0° 

3 

1,1, 3, 3, 1,1, 1,3, 3, 3, 5, 5, 5, 5, 5, 3, 5, 5, 7, 7, 5, 7, 1,7, 7, 1,1, 7  {28} 

o 

O 

oc 

4 

1,1, 1,1, 1,1, 1,3, 3, 3, 3, 3, 5, 5, 3, 5, 5, 5, 5, 5, 5, 5, 7, 7, 7, 7, 7, 1,1, 7  {30} 

0° 

5 

1,1, 1,3, 1,1, 1,1, 3, 3, 3, 5, 5, 5, 5, 5, 5, 5, 5, 7, 7, 7, 1,7  {24} 

0° 

6 

1,1, 1,1, 3, 1,3, 3, 1,3, 3, 3, 3, 5, 3, 3, 5, 3, 5, 5, 7, 5, 5, 5, 7, 5, 7, 5, 5, 7, 5, 7, 

1,7, 1,7, 7, 1,7, 1,1, 7  {42} 

0° 

7 

1,1, 1,3, 1,1, 1,3, 3, 3, 3, 3, 3, 3, 5, 3, 5, 5, 5, 5, 5, 7, 7, 7, 5, 5, 7, 7, 7, 7, 7, 1,1, 7 
{34} 

270° 

Results  of  a  similarity  measure  computing  for  analyzing  objects  using  variance  ratios  (5)  -  (7),  the  ranking 
coefficient  and  an  add-on  of  lacking  direction  vectors  weights  of  the  etalon  or  the  analyzing  object  by 
mathematical  expectation  meanings  are  located  in  the  Table  3.  It  should  be  stated  that  ranking  of  analyzing 
objects  (table  3)  in  decreasing  order  of  their  variance  ratios  has  the  same  distribution  for  utilized  variance 
ratios  (5)-(7). 


Table  3 


Computation  of  variance  ratios  value 


Object 

number 

rl  M 

72 

73 

1 

1,00000000 

1,00000000 

1,00000000 

2 

0,89260599 

0,91168540 

0,81373810 

3 

0,90017516 

0,91681997 

0,82529857 

4 

0,82261843 

0,86930268 

0,71510441 

5 

0,90868641 

0,92281044 

0,83854531 

6 

0,84103261 

0,87954314 

0,73972447 

7 

0,83581908 

0,87658396 

0,73266559 

X  Ar/i/n 

0,11415176 

0,08903634 

0,19070336 

It  follows  from  the  Table  3  that  variance  ratio  values  (7)  have  the  smallest  similarity  measures  for 
analyzing  objects  relative  to  values  computed  with  equations  (5)  and  (6).  Hence  the  variance  ratio  (7) 

^  Ar/3  /  n  =  0,19070336  is  more  sensitive  to  the  similarity  measure  of  analyzing  objects  relative  to 

distinguishing  [observable]  features  of  direction  vectors  weights  in  the  capacity  of  their  boundary 
representation. 

Boundary  representations  of  7  analyzing  objects  with  their  variance  ratios  are  depicted  on  the  Figure  5. 
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*%talon  -  100% 


4-71.510%: 


2  -  81.378%; 


C 


5  -  83.855%;  | _  _ J 

| - 1  7  -  73.267% 


6  -  73.972%; 


Fig.  5  Boundary  representations  of  analyzing  objects  with  their  variance  ratios 

Classification  results  of  analyzing  objects  by  threshold  values  of  their  variance  ratio  are  located  in  the 
table  4. 


Analyzing  objects  classification 


Table  4 


Object  numbers 

MCR,  % 

Object’s 

cluster 

number 

Cluster 

number 

Thresholds 
Similarity  measures 

i 

100,00 

i 

i 

99,953  -  100,000 

2 

81,373 

10 

2 

99,730  -  99,952 

3 

82,530 

10 

3 

99,400  -  99,729 

4 

71,510 

11 

4 

98,760  -  99,399 

5 

83,855 

9 

5 

97,560-98,759 

6 

73,972 

11 

6 

95,450-97,559 

7 

73,267 

11 

7 

92,810-95,449 

- 

- 

8 

89,040  -  92,809 

- 

- 

9 

83,850  -  89,039 

- 

- 

10 

76,990  -  83,849 

- 

- 

11 

68,270  -  76,989 

- 

- 

12 

57,630-68,269 

- 

- 

13 

45,750-57,629 

- 

- 

14 

31,080-45,749 

- 

- 

15 

15,850-31,079 

- 

- 

16 

0,800-  15,849 
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System  likeness  for  object  identifying  on  remote  sensing  images  is  shown  on  the  fig.  6. 


Fig.  6.  Object  identifying  system  likeness 


CONCLUSION 

The  sensitivity  of  variance  ratio  computing  is  increasing  up  to  5,74%  relative  to  utilized  variance  ratio  (5) 
by  forming  and  using  of  equinumerous  set  of  features  (adding  direction  vectors  weights  by  mathematical 
expectation  and  computing  the  ranking  coefficient  of  difference  between  etalon  vector  weight  and 
analyzing  object  vector  weight)  for  classification  of  analyzing  objects  on  remote  sensing  images  by  their 
boundary  representation  with  the  aid  of  variance  ratio  (7)  of  their  vectorial  models. 
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